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(57) Abstract 

Method of monitoring and/or diagnosing the operation- 
al state of a complicated, dynamic system constituted by a ma- 
chine or an industrial process installation, a plurality (M) of 
state variables being measured at a plurality (N) of instants for 
the assembly of a corresponding measurement data (XjjJ, 
which are processed statistically to form at least one refer- 
ence class corresponding to a specific operational state. 
There is used here a group of measurement data Qi\y) serv- 
ing as a learning group, which is adapted by linear projec- 
tion to a new coordinate system with fewer coordinate axes. 
In a corresponding way, current measurement or calculation 
data are adapted to the new coordinate system, whereupon 
the distance between the current measurement or calculation 
data, on the one hand, and the measurement data corre- 
sponding to the learning group, on the other hand, is deter- 
mined. If this distance exceeds a predetermined value, it is 
established which state variables that cause the deviation, by 
calculating the effect of the respective state variable on the 
variance along the respective coordinate axis in the new 
coordinate system. 
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A METHOD FOR MONITORING THE OPERATIONAL STATE OF A SYSTEM. 

The invention relates to a method of monitoring and/or 
diagnosing the current operational state of a complicated 
dynamic system constituted by a machine, apparatus or an 
industrial process installation, in accordance with the pre- 
amble to claim 1 . 

Such a method, intended for monitoring machine vibrations while 
utilizing frequency analysis as well as sampling and pattern 
recognition techniques is known from EP-A-84902732 . 1 , the 
reference class being frequency spectra and the pattern recog- 
nition and detection means being adapted such as to calculate 
the probability for each new frequency spectrum that it belongs 
to a class other than the reference class. An abnormal opera- 
tional state of the machine is determined when this probability 
exceeds a predetermined limit. 

According to the present invention, it is not only a question o 
analysing occuring vibrations, but a large number of measurable 
or calculable physical and possibly even logical parameters or 
state variables, which to a greater or less extent may be con- 
sidered as effecting the operational state of a complicated 
system . 



For example, in the operation and maintenance of complicated 
machines, e.g. jet engines, measurements of a large number of 
statistic and dynamic state variables are collected. 
Accordingly, being able to readily indicate an unfavourable 
combination of such variables is of the greatest importance 
here, not only with regard to limit values decided for each 
individual variable. The situation is similar for process 
installations, e.g. paper manufacture or energy conversion 
(such as in nuclear power stations). 
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There is further a desire of being able to calculate beforehand 
whether a proposed modification of a proven system will give 
acceptable results . 

5 A parmount object of the present invention is thus to monitor 

and/or diagnose, possibly by calculation performed beforehand, 
and with the aid of mathematical, statistical methods, the 
existing operational state of a complicated machine, apparatus 
or process installation and in such a case either identify a 
10 specific operational state (corresponding to a reference class) 

or determine a deviating operational state by simultaneously 
identifying what state variables cause the deviation* 

This object is achieved by the process of steps disclosed in the 
15 characterizing part of claim 1 . 

Statistic methods known per se are thus applied in accordance 
with the invention, these methods being used within so-called 
multivariate analysis, e.g. for rapidly discovering an undesired 
20 state in a machine and pointing out individual variables or 

groups of variables responsible for this state. The method may 
be applied either in real time, namely to incoming measurement 
values, or to data bases, where the measurement values have been 
stored . 

25 

The invention will know be explained in more detail below and 
with reference to the accompanying drawings, where 

Fig. A is a schematic diagram over the measuring and data equip- 
30 ment , its structure being known per se, which is used in the 

method according to the invention; 

Figs. 1-3 are flow diagrams corresponding to some application 
examples A, B and C ; 

35 
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Fig. la schematically illustrates a variant of the application 
example A; and 

Fig. 3a is a variant of the flow diagram in Fig. 3. 

5 

The invention is founded on the following procedures: 

The measurement values for different state variables Vi are 
obtained from suitable measuring sensors after amplification. 

10 filtering and digitizing Csee Fig. A). For rotating machines 

both static and dynamic state variables may be used, e.g. the 
dimensions, taken from a coordinate table, of the machine com- 
ponents, clearance and play between given parts, throw and im- 
balance measured in a balancing machine, as well as different 

15 types of vibrations. 

Measurement of M different state variables V ± is repeated at 

different instances N times. The results <X ik , i = 1, ...» H; 

k = 1, .... N), which form a matrix or table with the dimensions 

20 M x N can be presented in a N-dimens ional space as a cluster of 

N dots. The appearance of the cluster can be described, for 
example with the aid of a covariance matrix. If the dotB have 
been collected during two different states of the system they 
will form, with the greatest probability, two clusters in the M 

25 space, provided that among the state variables there are pre- 

cisely those which differentiate the current states in the 
system. A possibility is that the major part of the measurement 
data, which are derived from one or more known states of the 
system (so-called learning groups), are combined with a plu- 

30 rality of current measurement data which are to be tested <a 

test group) . 

For easy handling of this type of multivariate measurement data 
it is proposed in accordance with the invention that a co- 
35 ordinate transformation is made. For such a transformation 

to be utilizable it must: 
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- reduce the number of coordinate directions; 

- maintain the variance of the dot clusters; and 

- separate different dot clusters in the best way. 

A transformation found to meet these conditions is the so-called 
principal component method, which has been previously used in 
chemical data processing and in qualitative analysis of food- 
stuff. In this connection reference is made to the article 
"Pattern Recognition: Finding and using regularities in multi- 
variate data" by S. Wold, C. Albano, W.J. Dunn, K. Esbensen, S. 
Hellberg, E. Johansson and M. Sjostrom in the publication "Food 
Research and Data Analysis" Applied Science Publishers, London 
and New York, 1983. 



Applied to monitoring the state of machines, apparatus or 
industrial process installations, this principal component 
method means that multivariate measurement data X ik are 
approximated by linear projection to a new coordinate system 
with fewer coordinate axes (L; L < M) according to the 
20 relationship: 



L 

X iJc ■ Ai + Z B ia T ak + R ik 
a = 1 

where the new variables A, B and T describe the systematic part 
of the measurement data X ik , and where the residue matrix R ik 
describes the random part of X ik . It should be pointed out that 
the transformation separates measurement data into parameters 
(A, B) describing the state variables, on the one hand, and 
parameters (T) describing the system state at the respective 
measurement (k), on the other hand. The parameters T ak describe 
the location of the k-state in the M-dimensional dot cluster and 
are therefore used for discriminating state classes. 

Adaptation of the principal components takes place in conse- 
cutive steps, one component at a time. 
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The procedure starts by normalizing the measurement data matrix 
Xifc with the standard deviation S£ of the respective state 
variable V± , so that all variables will be of appoximately the 
same order of magnitude, the mean value Ai being substracted, 
5 i.e. the matrix X±^ - is formed. At this stage the first 

pr inc ipal component 

PCI = B ±1 T ik 

10 is adapted while using the so-called NIPAL's algorithm (Appendix 

A in the above mentioned article by Word et al>, which in short 
means that beginning with an attempt or qualified guess with 
respect to the vector B, the vector T is calculated and the 
convergence is checked, which is repeated until the desired 

15 accuracy has been achieved. The first residue matrix Rifci a 

x ik~ A i~ PC1 is then formed, which is used as the input matrix for 
determining the next (second) principal component, etc. 

The number of components (L) is determined by the character of 
20 the data matrix. The process of adaptation can be terminated 

when no systematic parts remain in the residue matrix R^fe. For 
practical reasons the process can also be terminated after PC2 . 
The state can then be depicted on a plane with coordinate axes 
corresponding to PCI and PC2 . A three-dimensional depiction can 
25 be constructed if the adaption is terminated after PC3 . Three- 

-dimensional depiction of measurement data are generally pre- 
ferred . 

It is of the greatest interest to find out how much different 
30 state variables contribute to the variance along each principal 

coordinate. In order to separate different state classes it is 
required, as mentioned above, that measurement data really do 
contain the state variables which are relevant for describing 
the classes. To decide this, the so-called modelling capacity 
35 P^ a of the state variable along the respective coordinate 

axis a (a = 1 , 2 , . . . , L) is defined: 
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?ia - 1 - SiaCRika> / "i <*ik - Ai> 

where 6i a is the standard deviation for the residue matrix Rifca 
5 and s^ is the standard deviation for the matrix CXifc - Ai). Pi a 

varies between 0 and 1. If P± a is close to 0, the variable 
does not affect the variance along the principal coordinate in 
question, and in principle can,, be ignored. The modelling capa- 
city can also be expressed as a percentage. In addition, the 
10 total modelling capacity for PCI, PC2 and PC3 is defined 

according to the following; 



'it = \/<Fil + A.2 + Pi3>/3 



15 Variables with a total modelling capacity Pjt close to 1 are 

denoted "strong" variables. Strong variables are relevant for 
identifying different states, particularly fault states. Know- 
ledge of these variables helps the operator to identify faults 
in machines, for example, or may be used in coaction with a 

20 knowledge base for automatic fault identification. 



Some practical examples A, B and C of the method in accordance 
with the invention are given below. 



25 Example A. Analysis of a data base with motor parameters 

( se Fig . 1 ) 

In connection with the inspection of a jet engine, about a 
hundred different static and dynamic variables in the engine are 
measured. Thus, e.g. tolerance and throw in different parts of 
30 the rotor-stator system are measured accurately. The measure- 

ments are placed in a data base. When an engine has been 
approved after inspection, these measurement values are stored 
together with the engine documents. On the other hand, when an 
engine has been rejected, attempts have been made in different 
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ways to search among about a hundred measurement variables for 
the one or ones which have caused the rejection. This is an 
insolvable problem using elementary stastics as the sole aid. 
The use of conventional packages for advanced statistic analysis 
5 is unfortunately too time-consuming for enabling their opera- 

tional use. 

However, in applying the inventive method, about ten data base 
posts from approved engines are used as a learning group. The 

10 post for the tested, rejected engine is added to the table 

(matrix). After principal coordinate transforming, the variables 
are ordered in rank with a modelling capacity according to size. 
In this way the one or ones of the state variables which are 
responsible with the greatest probability for rejection of the 

15 motor can be indicated. In combination with a knowledge base <a 

rule system) directions as to required technical measures can be 
worked out. If the method is applied before the final check in a 
test bench economies can be made in many cases, since the number 
of rejections after checking can be reduced. 

20 

A variation of this example is that the data vectors (the 
columns in the total data matrix) for the approved engines are 
divided into sections corresponding to different work stages J 
(see Fig. la) during assembly of the engine (j = 1, ....P). In a 

25 corresponding way the data vector is divided for the tested 

engine and the analysis is made of each terminated work stage of 
a growing data matrix. This gives the possibility of breaking 
off the costly assembly process as soon as it is found that the 
tested engine will belong to the group of non-approved engines . 

30 At this stage the problem can be solved by using another selec- 

tion of parts or further adjustment. 
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Example B. Checking blade distribution in rotor discs 
(see Fig. 2) 

Hitherto, the blades in a rotor disc are distributed according to 
simple rules decided by the manufacturer, with respect to the 
momentum weight of each" blade. Computer programs calculating the 
distribution of blades around the rotor disc have also begun to 
come into use. Such a program distributes the blades with 
reference to their momentum weights so that inbalance is 
minimized. The method in accordance with the present invention 
enables the checking of a blade distribution arranged according 
to experience or calculation by comparison with previous, 
successful blade distributions, preferably from the same engine. 
Information in a data base containing such distributions of 
momentum weights which have given satisfactory dynamic 
properties are compared with a proposed distribution. It iB here 
possible to obtain a measure of the quality of goodness of the 
distribution (distance to the group "successful distributions"), 
and an account of the blades causing possible negative dynamic 
properties in the rotor stage (strong variables) can thus be 
obtained. It is also possible to test several proposed distri- 
butions simultaneously . 

Example C. Monitoring and ad vising with the engine running on a 
test bench (see Fig. 3). 

When running an engine on a test bench, a large number M (usual- 
ly about a hundred) operational variables are measured and 
stored. In modern test bench equipments the positions of all 
manual controls handled by the operator are also read off. 
Hitherto, the automatic monitoring is arranged such that maximum 
and minimum limits are given for certain essential variables, 
and that the system gives an alarm for certain programmed com- 
binations of discrete and logical variables. The system thus 
gives an alarm for previously known undesired situations. 
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The inventive method can here be applied according to two 
different variations for an engine running on a test bench: 

Variation 1 . A large learning group is formed containing 
5 2-3 measurements for each permitted operational situation 

(corresponding to a specific reference class). Each new measure- 
ment ("present state" in Fig. 3) is continuously added to the 
table or data matrix thus obtained. The new measurement (a dot 
in the M-diment ional space) is tested with respect to its 

10 distance to the learning group. If the new dot is outside those 

of the learning group, a deviating stage is established and 
possible strong variables are accounted for. If measures taken 
by the operator give rise to strong variables the system can, 
with support of a knowledge base, give the operator advice and 

15 information about suitable measures in this situation. Possible 

operator errors can be discovered and dealt with in this way. 

Variation 2 . Several learning groups are formed, one for each 
interesting operational situation. Each learning group contains 

20 about ten measuring points coming from this situation. Certain 

undesired situations, e.g. faults or errors, can be stored among 
the learning groups. Each new measuring point, or several con- 
secutive incoming measuring points can then be tested against 
the learning groups. By searching for the best adaptation to a 

25 learning group, the operational situation can be identified with 

great certainty. A simultaneous account of strong variables can, 
in combination with a knowledge base, be used for achieving a 
diagnosis and for giving advice with respect to possible 
measures . 

30 

Both variations can be used both in real time as well as 
' subsequently for analysis of collected measurement data. Both 
variations can further be used on all types of machines, 
apparatus or industrial process installations, where a large 
35 number of operational variables make the situation difficult for 

a human being to comprehend. 
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This method can be made even more effective (see Fig. 3a) 
if the data vector in question ("present state") is transferred 
to an autoassociative neural network, e.g. of the type ART 
(Adaptive Resonance Theory Network), as decribed in Carpenter, 
G.A., Grossberg, S., M A Massively Parallel Architecture for a 
Self -Organizing Neural Pattern Recognition Machine", Computer 
Vision, Graphics and Image Processing 37, pp 54-115, 1987. The 
network is trained with typical data from normal operational 
conditions. The remaining part of the data flow is only acti- 
vated on the discovery of a new category, for then carrying out 
the complete statistical analysis. In this way the statistical 
analysis does not need to be made during normal operational 
conditions . 
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CLAIMS 

1. Method of monitoring and/or diagnosing the current opera- 
tional state of a complicated, dynamic system constituted by 
5 a machine with a plurality of components, an apparatus or an 

industrial process installation with a plurality of state 
variables specific to the process, such as temperatures, flows, 
etc, a plurality (M) of physical and possibly also logical state 
variables (V ± ,i = 1,2, . M> for said system being measured at 

10 a number (N> of instances <t k , k = l, 2, . N) for making up a 

corresponding number (MxN) of measurement data <Xik; i = 1, 2, 

M; k = 1,2, N), which are processed according to 

statistical methods in a computer for forming a plurality of 
reference classes corresponding to specific operational states 
15 of the system, the current measurement or calculation data <Xiki 

i = 1,2, M) being compared with said reference class or 

classes for diagnosing the corresponding, current operational 
state of the system, characterized in that 
- a group of measurement data (XikJ * = 1*2, ...,M; 
20 k = 1,2, N';N' < N) derived from the measurement of 

said state variables serves as a learning group for forming 
the respective reference class, these measurement data <X^k> 
being adapted in a manner known per se by linear projection 
to a new coordinate system with fewer coordinate axes 
25 <L;L < M) while retaining the variance of the measurement 

data and the selectivity of different operational states, 
according to the relationship 
L 

x ik = A ± + Z B ia T ak + Riki i -l,2 f ...,M;k =1,2...,N' 
30 a=l 

where Aj_ is the mean value of the respective state variable 
<Vi> ; 

Bi a and T a k are a reduced number of parameters in the new 
35 coordinate system essentially corresponding to the 

systematic part of the measurement data; 
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R ik is a residue matrix essentially corresponding to the 
random part of the measurement data; 
- said current measurement or calculation data (X^; i = 

= 1,2, ... M; N'<k<N) are adapted in a corresponding manner 
to said new coordinate system, 

the current operational state of the system, possibly after 
checking in an autoassociat ive network as to whether the 
current operational state belongs to any reference class, is 
compared with said specific operational state by determining 
the mutual distance between said current measurement or 
calculation data, on the one hand, and the measurement data 
corresponding to the respective learning group, on the other 
hand, in the new coordinate system, and 

in the case where said possible check has not been carried 
out and said distance falls below a predetermined value, the 
corresponding specific operational state of the system is 
identified and, in the case where said distance exceeds a 
predetermined value, it is determined what state variables 
and, on the basis thereof, what component(s) or state 
variables(s) specific to an industrial process^ respectively, 
that cause the deviation by calculating the effect of the 
respective state variable on the variance along the 
respective coordinate axis in the new coordinate system. 

2. Method as claimed in claim 1, characterized in 
that said adaptation of measurement data CX^) to the new co- 
ordinate system is carried out in steps, the respective co- 
ordinate axis in the each step (a; a = 1,2, . L) being 
selected such that the variance for said parameters <Bi a ,T a j c ) 
along this coordinate axis is minimized. 

3. Method as claimed in claim 2, characterized 
in that said stepwise adaptation is carried out according to 
NIPAL's algorithm. 
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4. Method as claimed in claim 2 or 3, characteriz- 
ed in that the new coordinate system has one, two or three 
coordinate axes (L=l,2 or 3). 

5 5. Method as claimed in any one of claims 2-4, charac- 

terized in that said calculation of the effect on the 
variance along the respective coordinate axis caused by the 
respective state variable is carried out on the basis of the 
so-called modelling capacity CPi a J a = 1, . . L) of the 
10 respective state variable according to the formula: 

Pia = 1 " s ia CR ika )/ Si (X ik - Ai>; 

where s^ a and Sj are the standard deviations of Rika 

and Ai , respectively; and 

R ika is the residue matrix after the respective 

15 adaptation step (a); 

this modelling capacity (Pi a > being able to assume values, 
between 0 and 1 and the total modelling capacity <Pit> of tne 
state variable being calculated as a combination of the values 
of the different adaptation steps <a> so that a total modelling 

20 capacity between 0.5 and 1 points out the state variable in 

question as contributing to said deviation. 

6. Method as claimed in claim 5, characterized in 
that the total modelling capacity <Pit^ of the respective state 
25 variable is calculated according to the formula: 



2 2 
P ± l + ... + P iL 



'it 



30 



7. Method as claimed in any one of the preceding claims, 
characterized in that said measurement or cal- 
culation data (Xifc) are normalized with the standard deviation 
35 (s^) of the respective state variable, calculated on said 

learning group, before said adaptation to the- new coordinate 
system . 
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8. Method as claimed in any one of the preceding claims for 
analysing previously made measurement data stored in a data base 
and relating to several systems of the same type, said systems 
having been approved of or rejected on inspection or test 
5 operation, characterized in that a plurality of 

measurement data posts derived from approved systems are used 
as a learning group, and one or more measurement data posts 
derived from rejected systems are used as said current measure- 
ment data, there being produced a list of the state variables 
10 which have, with a certain probability, contributed to rejection 

of the system. 



9, Method as claimed in any one of claims 1-8 for monitoring 
different operational states in test running of complicated 

15 systems, which are at least partially controlled manually by an 

operator, characterized in that a plurality of 
measurement data posts derived from permitted operational states 
are used as one or more learning groups, which are successively 
supplemented with current measurement data derived from the test 

20 run, deviating operational states and state variables causing 

these being accounted for continuously for diagnosis and 
possible measures to be taken by the operator. 

10. Method as claimed in any one of claims 1-7 for calculating 
25 beforehand whether a proposed modification to a tested system 

will give acceptable results, characterized in 
that measurement data posts, derived from a plurality of tested 
and approved systems of a similar kind, are used as a learning 
group, whereas computer calculated data corresponding to said 
30 proposed modification are used as said current data. 
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